Low-Mass Dileptons at the CERN-SpS: Evidence for Chiral Restoration? 
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Using a rather complete description of the in-medium p spectral function - being constrained 
by various independent experimental information - we calculate pertinent dilepton production rates 
from hot and dense hadronic matter. The strong broadening of the p resonance entails a reminiscence 
to perturbative qq annihilation rates in the vicinity of the phase boundary. The application to 
dilepton observables in Pb(158AGeV)+Au collisions - incorporating recent information on the hadro- 
chemical composition at CERN-SpS energies - essentially supports the broadening scenario. Possible 
implications for the nature of chiral symmetry restoration are outlined. 
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I. INTRODUCTION 

In the low-energy sector (q 2 < 1 — 2 GeV 2 ) the proper- 
ties of Quantum Chromodynamics are governed by (ap- 
proximate) chiral symmetry and its dynamical breaking. 
In hot and dense matter as created through energetic col- 
lisions of heavy nuclei chiral symmetry is believed to be 
restored, associated with the disappearance of the chi- 
ral condensate and a substantial reshaping of the low- 
lying hadronic spectrum. As the chiral condensate does 
not represent a physical observable, the detection of sig- 
natures for chiral restoration has to focus on medium 
effects in hadronic properties. Dileptons as penetrat- 
ing probes provide the opportunity to study rather di- 
rectly the in-medium modifications of light vector mesons 
through their decays V — > l + l~ . 

Measurements of low-mass dilepton spectra in heavy- 
ion collisions at CERN-SpS energies [^]J^] have revealed 
appreciable radiation beyond final-state hadron decays. 
However, the spectral shape of this additional contri- 
bution - stemming from the interaction phase of the 
hadronic fireball - is not easily accounted for. In par- 
ticular, using the vacuum vector meson line shapes, the 
enhancement around Mu ~ (0.2 — 0.6) GeV remains un- 
explained. At the same time the yield in the p/u; mass 
region tends to be overestimated. The central question 
then is whether one has possibly observed signatures of 
(the onset of) chiral symmetry restoration, and, if so, 
how it manifests itself in the vector channel. 

Several theoretical approaches have been pursued to 
study medium effects in dilepton production rates from 
hot/dense matter, e.g., the 'dropping rho mass' scenario 
(within a mean-field type treatment) [[|-||, the chiral re- 
duction formalism (based on chiral Ward identities in 
connection with low-density expansions) chiral La- 
grangian frameworks or many-body-type calculations 
of in- medium vector- meson spectral functions p-pT| . As 
a further complication, the quantitative impact of in- 
medium effects on the final spectra depends on another 
important ingredient, namely the modeling of the space- 
time evolution of the global heavy-ion collision dynam- 
ics. Here, transport calculations seem to give substan- 
tially larger dilepton yields than hydrodynamic simula- 



tions. A reason for this discrepancy might be related 
to an overpopulation of pion phase space (enhancing the 
7T7T — » p — > 11 annihilation contribution). 

In the present article we will focus on the in-medium 
spectral function approach for the elementary dilepton 
production processes, presenting its most recent devel- 
opments, thereby noting a surprising reminiscence to 
(perturbative) quark-based calculations. We then ap- 
ply the medium-modified dilepton rates to the calcula- 
tion of spectra from iti: —> p — > ee annihilation in 30% 
central Pb(158 AGeV)+Au collisions and perform a de- 
tailed comparison to various projections of the recent 
data from the CERES /NA45 collaboration. Special care 
is taken to adopt a realistic temperature/density profile 
of the collision dynamics using a thermal fireball expan- 
sion with a particle composition consistent with recent 
hadro-chemical analyses at SpS energies The 
subsequent evolution towards thermal freezeout is deter- 
mined by assuming entropy as well as pion- and kaon- 
number conservation which induces the build-up of finite 
pion- /kaon-chemical potentials . 



II. PRODUCTION RATES IN HOT/DENSE 
MATTER 

The radiation of dileptons from a hot and dense 
medium characterized by a temperature T and baryon 
chemical potential /j,b, 
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is governed by the thermal expectation value of the (re- 
tarded) electromagnetic current-current correlator, n em 
(f B : Bose function). The latter can be written in terms 
of its longitudinal and transverse projections as 

U em {q ,q) = -[U^ m {q ,q)+2Ul m {q ,q)] , (2) 

both depending separately on energy and momentum. 
Below momentum transfers of about 1 GeV, the hadronic 
part of the e.m. current can be accurately saturated by 
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the light vector mesons within the well-established vec- 
tor dominance model (VDM). The correlator is then ex- 
pressed through the imaginary parts of the vector meson 
propagators (= spectral functions) as 
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In the following we will neglect the small contributions 
from the isoscalar part and concentrate on the p meson 
which plays by far the dominant role for the time scales 
involved in heavy-ion reactions. 

Along the lines of our earlier analyses || the p propa- 
gator in hot hadronic matter, 
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is evaluated in terms of various contributions entering 
its in- medium selfenergy £ p . The meson gas effects 
(encoded in £ p m) are accounted for following Ref. 
through interactions with the most abundant thermal it, 
K and p mesons saturated with a rather complete set of s- 
channel resonances R up to 1.3 GeV. The interaction ver- 
tices have been constrained by both hadronic (R — > Pp) 
and radiative (R — > P-f) decay branchings. Also, the 
Bose enhancement for the in-medium pin: width has been 
included in S pT7r . 

The p modifications in nuclear matter are based on 
the model constructed in Refs. ]s|, |l1||l7t , where direct 
pN -> B interactions (B—N, A, V(1520), A(1700), 
V(1720), . • ., encoded in £ p b) as well as pion cloud 
modifications through irNN -1 and 7rAV _1 excitations 
(entering £ p7r7r ) were calculated. Again, the model pa- 
rameters have been thoroughly constrained by analyz- 
ing photoabsorption spectra on nucleons and nuclei as 
well as ttN — > pN scattering data. The latter enforce a 
rather soft form factor on the irNN vertex. Also, since 
the VDM is less accurate in the baryonic sector, an im- 
proved VDM coupling has been employed for the 
(transverse) jN — > B transition form factors so that 
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with 
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cf. Refs. |j,[l(| (here tb — 0.7, the ratio between the 
actual jNB coupling and its value in the simple VDM). 
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FIG. 1. Real (lower panel) and imaginary (upper panel) 
part of the p propagator in hot and dense hadronic mat- 
ter at fixed baryon chemical potential \xb = 330 MeV and 
temperatures T = 180 MeV (corresponding to nucleon and 
total baryon densities of qn = 0.67go and qb = 2.6£o), 
T = 150 MeV (qn = 0.25^ o , Qb = 0.72£> o ) and T = 120 MeV 
(pjv = 0.06p , Pb = 0.13p ). 

At finite temperatures two additional features appear 
in the baryonic sector. Firstly, the Fermi-Dirac distribu- 
tion functions are replaced by thermal ones Jn||| 1 which, 
at the temperatures of interest here (T ~ 150 MeV), 
are substantially smeared as compared to the zero- 
temperature 0-functions. Secondly, a large fraction of 
the nucleons is excited into baryonic resonances which, in 
turn, are no longer active in nucleon-driven effects. How- 
ever, as conjectured in Ref. pQ] , a baryonic resonance 
B\ might still have excitations of type pi^-Bf 1 built on 
it. The particle data table |21| indeed supports this hy- 
pothesis: e.g., the £(1670) (which is a well-established 
four-star resonance with spin-isospin IJ P — l| ), when 
interpreted as a p£ (or pA) 'resonance', very much re- 
sembles the quantum numbers and excitation energy 
(AE ~ 500 - 700 MeV) of the pN -> AT(1520) transi- 
tion. In addition, the branching ratio of £(1670) decays 
into 'simple' final states such as NK, T,ir or A7r is sig- 
nificantly less than 100% . Similar excitations on non- 
strange baryonic resonances are more difficult to identify 
as the latter themselves decay strongly via pion emis- 
sion (i.e., the Bi — > B\p decay is immediately followed 
by further Bi — ► Nix and p — > 7T7t decays). Neverthe- 
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less, from pure quantum numbers it is tempting to as- 
sociate, e.g., A(1930)A _1 or iV(2080)7V(1440)" 1 excita- 
tions with S'-wave 'Rhosobar' states. We have conser- 
vatively estimated their coupling constants to give rea- 
sonable branching fractions into pB\ and 7B1 (the latter 
are typically in the range of 0.2-0.7 MeV). Due to the 
large total widths involved, their impact on the p propa- 
gator is weaker as compared to excitations on nucleons at 
equal densities; e.g., at T=150 MeV and p,s = 436 MeV, 
where, in chemical equilibrium, about 2/3 of the baryons 
are in excited states, the additional broadening from the 
calculated -E^-Bi" 1 excitations as compared to BN^ 1 ones 
amounts to about 40%. Concerning the pion cloud mod- 
ifications, we do not explicitly compute the 'Pisobars' on 
excited resonances but approximate their effect by using 
an effective nucleon density g e ff — qn + 0.5qb* ■ 

The final result for the p spectral function in hot 
hadronic matter is displayed in Fig. |l| at fixed three- 
momentum and chemical potentials p^ = 0, ps = 
330 MeV for pions and baryons, respectively. One ob- 
serves a strong broadening with increasing temperature 
and density. Comparing to the pure meson gas results of 
Ref. |l5| (cf. Fig. 3 therein), it becomes clear that espe- 
cially the low-mass enhancement around M ~ 0.4 GeV 
is largely driven by baryonsf]. 

The pertinent three-momentum integrated dilepton 
production rates are shown in Fig. |^. At moderate tem- 
perature and density (upper panel), one still recognizes 
a remnant of the p peak, which, however, is entirely 
wiped out under conditions expected to be close to the 
phase boundary (lower panel). This provokes a compari- 
son with quark-gluon based rate calculations, which, for 
simplicity, have been evaluated in terms of lowest order 
0(aoP s ) qq annihilation |22j|, 
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with x± = exp[-(g ± q)/2T], y = exp[-(q + Mg)/ T ]- 
In the vacuum the qq rates are known to coincide 
with the hadronic description for invariant masses M > 
1.5 GeV as marked by the cross section ratio cx(e + e~ — » 
hadrons) /a(e + e~ — ► p + p) for the inverse process of 
e + e~ annihilation (up to additional resonance struc- 
tures in the vicinity of heavy-quark thresholds such as 
cc). It seems that the in-medium hadronic rates in- 
deed approach the partonic ones rather quickly leading 



to an approximate agreement at the highest tempera- 
tures/densities for masses of about 0.5-1 GeV (the de- 
viation at masses M > 1 GeV is due to the incomplete 
description of the vector correlator in vacuum which docs 
not include more than two-pion states; the discrepancy 
at low masses might be reduced once higher order as cor- 
rections are included, in particular soft Bremsstrahlung- 
type processes) . A tempting interpretation of this behav- 
ior is a lowering of the 'quark-hadron duality threshold' 
in hot and dense hadronic matter. We will return to this 
issue in the discussion in Sect. 
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FIG. 2. Comparison of hadronic and partonic dilepton pro- 
duction rates (integrated over 3-momentum) in hot and dense 
matter at fixed baryon chemical potential (and equivalent 
quark chemical potential p, q =p,B/3) for T=150 MeV (upper 
panel) and T=180 MeV (lower panel); dashed and full lines: 
tttv annihilation using the free and the in-medium p spec- 
tral function, respectively; dashed-dotted lines: lowest-order 
quark-antiquark annihilation. 



III. LOW-MASS DILEPTON OBSERVABLES IN 
PB(158 AGEV)+AU 



*With respect to the conditions realized in URHIC's at the 
full SpS energy (160-200 AGeV) this statement has to be 
taken with care as finite pion chemical potentials might arise 
towards freezeout which enhances the contributions from the 
meson gas. 



For a sensible application of the in-medium vector cor- 
relator to calculate dilepton production in URHIC's a 
realistic temperature and density profile is required. Mi- 
croscopic transport or hydrodynamical simulations, e.g., 
have proven very successful in describing the final hadron 
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observables. However, the associated dilepton yields, 
based on identical production rates, may differ appre- 
ciably H|, cf., e.g., Refs. [|§ and §|. As we will show 
below, a possible origin of this discrepancy might be re- 
lated to the build-up of a finite chemical potential, which 
is usually not included in most hydrodynamical calcula- 
tions to date (see, however, Ref. pE - 



Recent hadro-chemical analysis |12[Jf 3f| of a large body 
of hadronic heavy-ion data has shown that the finally 
observed particle abundances at SpS energies are con- 
sistent with a common chemical freezeout at temper- 
atures/baryon chemical potentials around (T,fis)ch — 
(175,270) MeV. In the subsequent expansion and cool- 
ing the system is still strongly interacting with elas- 
tic collisions maintaining thermal equilibrium until the 
thermal freezeout. The absence of pion-number chang- 
ing (inelastic) reactions (which is well supported phe- 
nomenologially from the inelasticities in ttt: scattering) 
then induces a finite pion chemical potential [ Q| . To 
incorporate these features, we here take recourse to a 
simple expanding thermal fireball model. We start by 
generating a (non-interacting) resonance gas equation of 
state of hot hadronic matter including the lowest 32 (16) 
mesonic (baryonic) states. Imposing entropy and baryon 
number conservation determines the temperature depen- 
dence of the baryon chemical potential, pb(T). The addi- 
tional requirement of pion- [kaon-] number conservation 
entails a finite p-n{T) \pk{T) ~ p^(T)], with the vari- 
ous baryon chemical potentials kept in relative equilib- 
rium (with respect to strong interactions), e.g., p\{T) = 
Pn(T) + p^fT), etc.. With an entropy per baryon of 
s/qb = 26, our trajectory is found to pass through 
(T,fx N ) = (175,250) MeV, compatible with the exper- 
imentally deduced point |12[ (here, a small fi w ~ 20 MeV 
is needed to be consistent with the final pion-to-baryon 
ratio of about 5:1 at SpS energies). Towards thermal 
freezeout at (T,p N ) fo ~ (110 - 120,415 - 450) MeV, 
p-niT) and pk{T) increase approximately linearly to 
around 80 and 110 MeV, respectively. Finally we need to 
introduce a time scale to obtain the volume expansion. 
We approximate the latter by a cylindrical geometry as 

V${t) = 2 (z + v z t + \a,t 2 ) 7T (ro + \a±t 2 ) 2 (8) 

employing two firecylinders expanding in the ±z direc- 
tion. Guided by hydrodynamical simulations |26[ the 
primordial longitudinal motion for Pb(158 AGeVj+Au 
reactions is taken to be v z — 0.5c, and the longitudinal 
and transverse acceleration are fixed to give final veloc- 
ities v z (tf ) ~ 0.75c, v±(tf ) ~ 0.55c as borne out from 
experiment (this, in turn, requires fireball lifetimes of 
about tf = 10 — 12 fm/c and implies transverse expan- 
sion by 3-4 fm, consistent with HBT analyses |27]]). The 
parameter ro denotes the initial nuclear overlap radius, 
e.g., ro = 4.6 fm for collisions with impact parameter 
b = 5 fm and Nb — 260 participant baryons. The pa- 
rameter zo is equivalent to a formation time and fixes the 
starting point of the trajectory in the (T, /ijy) plane. 



With the such specified space-time evolution for ^30% 
central Pb(158 AGeV)+Au collisions the dilepton yields 
from in-medium radiation are obtained as 



d 2 N ee 
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xF(M,q;p B {t),T(t)) Acc{M,q) 



(9) 



with Acc(M, q) accounting for the experimental accep- 
tance cuts (as well as an approximate mass resolution) in 
the '96 CERES data. The backward firecylinder is placed 
at a rapidity y — 2.6 to ensure the correct charged par- 
ticle multiplicity as quoted by CERES. The electromag- 
netic 'transition form factor' T contains the in- medium 
p spectral function as described above, where the addi- 
tional effects of the finite meson chemical potentials have 
been implemented in Boltzmann approximation. For ini- 
tial conditions (T, £>s)im = (190 MeV,2.55go) and a freeze- 
out at (T, g B ) / =(115 MeV,0.33g ) the fireball lifetime 
amounts to about i/ =ll fm/c, see also above (reducing, 
e.g., (T, qb)%ti% to (180 MeV,1.92p ) reduces the dilepton 
signal from the fireball by about 15%). 
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FIG. 3. 1996 CERES/NA45 inclusive dielectron invariant 
mass spectra from 30% central Pb(158 AGeV)+Au colli- 
sions [E8| compared to the hadronic cocktail [0 (without p 
decays, dashed-dotted line) and to the cocktail plus tvtt annihi- 
lation from the hadronic fireball using either the free (dashed 
line) or the in-medium (solid line) p spectral function from 
Eq. (§). 

In addition to the in-medium radiation described by 
Eq. (||) the finally observed dilepton spectra contain 
a sizeable contribution from free hadron decays after 
freezeout, the so-called 'cocktail'. For this we employ the 
most recent evaluation of the CERES collaboration [^9| , 
which is generated from particle abundances based on 
the chemical freezeout of Ref. |l2j being consistent with 
our (T,pb) trajectory. As in our previous works || the 
p-meson decays have been excluded from the cocktail as, 
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owing to the short lifetime of the p, they are accounted 
for through the in- medium part, Eq. (^), in the vicinity 
of the (idealized) 'freezeout'. 

The final results for the inclusive invariant mass spec- 
tra, shown in Fig. ||, demonstrate that the use of the in- 
medium spectral function leads to reasonable agreement 
with the '96 CERES data, which cannot be described as- 
suming vacuum p properties. This is in line with the con- 
clusions of our earlier analysis || , but we emphasize again 
that the results presented here are based on a much im- 
proved understanding of both the microscopic p spectral 
function (being constrained by a large body of indepen- 
dent data) as well as the space-time evolution dynamics 
(consistent with chemical freezeout analyses, without the 
use of an overall normalization factor, etc.). As compared 
to Ref. H the smaller pb in the early stages, the finite 
p n in the later stages, and the more complete assessment 
of the meson gas effects jl5| reduces the dominant role 
of the baryonic contributions to about equal importance 
as the mesonic ones. Another noteworthy point is that 
the overall dilepton yield exceeds what has been found in 
the hydrodynamical calculations of, e.g., Ref. (2j| (where 
/iyr = has been used throughout) by about a factor of 
~2.5, thereby essentially resolving the above mentioned 
discrepancy to transport calculations ||||. 
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FIG. 4. Same as Fig. [| but split into two transverse mo- 
mentum bins |2j| below (upper panel) and above (lower panel) 
q t =0.5 GeV. 



A more detailed insight into the nature of the low-mass 
enhancement is revealed by inspecting its dependence on 
the transverse pair momentum qt- Fig. |] shows the in- 
variant mass spectrum seperated into two transverse mo- 
mentum bins below and above qt — 0.5 GeV. Clearly, the 
excess of the data over the free calculation is most pro- 
nounced in the small momentum bin which is properly 
reflected in our in-medium calculations. The same be- 
havior is exhibited in another projection of the data in 
Fig. ^, where transverse momentum spectra in various 
invariant mass bins are displayed. 





| M ee <0.25GeV; 

At 


hK, 0.25 

-\.-' v s 


<M ee <o.e 


8GeV^ 


CE 
J>_». 


I I I I I I | I I I I I I I I I ; 

1RES/NA45 Pb+Au ; 
N 0.68<M aB <0.88GeV ; 

» "^54 - 


< 




50 





.88Ge 


V<M e8 ! 



0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 

q, ee [GeV] q, ee [GeV] 



FIG. 5. 1996 CERES/NA45 transverse momentum spec- 
tra for e + e~ pairs in 30% central Pb(158 AGeV)+Au colli- 
sions Line identification as in Fig. []. 

IV. SUMMARY AND DISCUSSION 

Based on realistic hadronic interaction vertices we have 
shown that the in-medium p spectral function, evalu- 
ated within standard many-body techniques, undergoes 
a strong broadening in hot and dense matter. At compa- 
rable densities the baryonic effects prevail over mesonic 
ones due to the stronger nature of the meson-nucleon 
interactions. Corresponding thermal dilepton produc- 
tion rates exhibit a remarkable tendency to approach the 
lowest-order perturbative qq annihilation with increas- 
ing temperature and density. This might indicate the 
emergence of 'quark-hadron duality' towards the phase 
boundary down to rather low mass scales of 0.5 GeV, 
made possible through large imaginary parts resummcd 
in the Dyson equation of p propagator. The medium 
effects might thus be interpreted as a lowering of the 
'duality threshold', which in free space is located around 
1.5 GeV. An unambiguous identification of chiral restora- 
tion, however, requires the simultaneous treatment of the 
axial correlator (ai channel) which has to merge with 
the vector correlator which we have focused on here. Al- 
though the perturbative qq rate implies restored chiral 
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symmetry, one has to keep in mind that (thermal) non- 
perturbative effects below 1 GeV are not yet well under 
control. The argument can be made more rigorous in 
the 1-1.5 GeV region, where the lowest-order in temper- 
ature mixing between vector and axialvector correlators 
suffices to establish a full degeneracy between them as 
well as the perturbatively calculated partonic result. 

We have furthermore demonstrated that our model for 
the in-medium p spectral function is in line with current 
low-mass dilepton measurements at the CERN-SpS. For 
that we have employed a simple fireball expansion model 
consistent with recent hadro-chemical analysis and im- 
posing effective pion-number conservation between chem- 
ical and thermal freezeout via finite pion chemical poten- 
tials (the latter seem to resolve a large part of the discrep- 
ancy in the total dilepton yield between transport and 
chemical equilibrium-based hydrodynamic simulations). 
The apparent depletion of the in-medium dilepton yield 
in the p/tu region, as well as the enhancement below can 
be accounted for. This also holds for the low-q t nature of 
the excess. Upcoming high resolution measurements as 
well as the more baryon-dominated 40 GeV run will put 
the model predictions under further scrutiny and can be 
expected to advance our understanding of chiral symme- 
try restoration in hot and dense hadronic matter. 
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